Twist-bend instability for toroidal DNA condensates 
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We propose that semiflexible polymers in poor solvent collapse in two stages. The first stage 
is the well known formation of a dense toroidal aggregate. However, if the solvent is sufficiently 
poor, the condensate will undergo a second structural transition to a twisted entangled state, in 
which individual filaments lower their bending energy by additionally orbiting around the mean 
path along which they wind. This "topological ripening" is consistent with known simulations and 
experimental results. It connects and rationalizes various experimental observations ranging from 
strong DNA entanglement in viral capsids to the unusually short pitch of the cholesteric phase of 
DNA in sperm-heads. We propose that topological ripening of DNA toroids could improve the 
efficiency and stability of gene delivery. 

PACS numbers: 64.70.Nd, 87.15.He, 61.30.Pq 
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Single polymers collapse from a random coil confor- 
mation to a dense state once the solvent gets sufficiently 
poor 0| . For a flexible chain the condition of minimal 
surface energy yields an approximately spherical glob- 
ule, but for semiflexible polymers the situation is more 
complex 0: The local structure of the dense phase then 
consists of essentially straight chains with a basically par- 
allel alignment, in order to minimize bending energy and 
maximize density, respectively. Such a state can be char- 
acterized by a smooth field of tangent vectors, but in 
the spherical case this field must have at least two en- 
ergetically unfavorable defects on the surface (you can't 
comb a sphere). However, for a torus many defect-free 
fields are possible. Indeed, DNA, the probably best stud- 
ied semiflexible polymer, readily forms beautiful nanotori 
after adding any one of a variety of possible condensing 
agents (like polyethylcnglycol (PEG), multivalent coun- 
terionSjOr bundling proteins) to a dilute solution of DNA 
chains 3]. These tori are surprisingly monodisperse, hav- 
ing a radius comparable to the persistence length of DNA 
(rs 50 nm) basically independent of the condensation 
method Q. 

Consider such a condensate, in which the chain is 
wound up like a garden hose to form a torus with ax- 
ial and tubular radii r\ and r2, respectively. Since r\ 
is the average radius of curvature of the chain, a sim- 
ple scaling analysis balancing a bending energy Ajr\ 
per unit length of polymer and a surface energy a per 
unit area of the torus yields Q t\ ~ (a/A)~ 2 / 5 V 1 / 5 and 
r 2 ~ (cr/A) 1 / 5 !/ 2 / 5 , where the chain volume V ~ T\t\ as 
well as the packing density are assumed constant. Hence, 
the aspect ratio £ = ri/ra ~ (cr/A) -3 / 5 !^ 1 / 5 shrinks 
if a or V increase (i e., if the solvent gets poorer or 
the chain longer), and thus the torus "fattens". In this 
case it is no longer justified to calculate the bending en- 
ergy with some average radius of curvature (p) = r\. In 
fact, since (p~ 2 ) > (p)~ 2 (by virtue of Jensen's inequal- 
ity), the actual curvature energy should be larger. How- 
ever, the same argument indicates that the condensate 



can lower its bending energy by redistributing curvature 
more evenly. In this letter we demonstrate that indeed 
below a critical aspect ratio £ c (or above a critical surface 
tension a c ) the system spontaneously relaxes bending en- 
ergy by twisting the bundle of polymer strands. 

Indirect indications of such a twisted state can be 
found in computer simulations Q. Analyzing Cryo- 
EM experiments on DNA toroids [f| Hud et al. pro- 
posed non-local DNA winding and equidistribution of 
bending in toroids 0. Besides indications for the non- 
trivial local organization of DNA in toroids there are 
several lines of evidence for a global, topological, non- 
triviality coming from in vivo studies: Certain bacterio- 
phages, whose DNA is (due to a genetic modification) no 
longer attached to their nucleocapsid, display unusually 
strong knotting of the genome pfl. And the chirality of 
the highly confined sperm-chromatin is surprisingly pro- 
nounced, with a pitch 10 times shorter than in vitro 
These findings make us wonder whether there is a con- 
nection to the topological ripening we will now discuss. 

Let us begin our quantitative analysis of the situation 
by neglecting the connectivity of the chain. More specifi- 
cally, we will first formulate a local theory which is based 
on the above mentioned nematic field n of unit tangent 
vectors The path of the actual polymer will later 
be recovered as an integral curve of this flow field, and 
its global topological properties can then be studied. The 
elastic energy e per unit volume, describing the deviation 
from perfectly parallel alignment, is the Frank-Oseen free 
energy of a uniaxial nematic liquid crystal 10] 
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+ ^( n .(Vxn)) 2 + ^(nx(Vxn)) 2 , 

(1) 

where the three terms correspond to splay, twist, and 
bend deformations, respectively. Assuming the conden- 
sate to behave like an incompressible liquid the first 
term (splay) must vanish identically in order to maintain 
a constant polymer density throughout the condensate 
[llj . while the other two terms divide the elastic energy 



FIG. 1: Illustration of the flow field on a toroidal condensate 
which features additional twist. The aspect ratio is £ = 1.5. 

between themselves. The total energy is of course the 
integral of Eqn. over the torus volume. 

It is convenient to discuss this situation in suitable 
toroidal coordinates {r, ip}, defined by 

x = (n —r cos i?) cos <p , (2a) 
y = (n — r cos §) sin ip , (2b) 
z = rsmi3 . (2c) 

The nematic field is now represented as n = n r e r + 
n$etf + n v e v (where the &i — dtr/\dir\ are the toroidal 
unit tangent vectors). When inserting this into Eqn. 
and integrating over the volume, we obtain the elas- 
tic energy as a functional of n r (r 7 d,ip), n$(r, -d, ip), and 
rij,(r, ip). At the surface of the torus the boundary 
condition n r (r = r-i) = must hold, and owing to rota- 
tional symmetry we will henceforth make the (nontrivial 
but very reasonable) assumption that none of the coor- 
dinate functions depends on ip. These steps reduce the 
task of finding the optimal polymer winding to a two- 
dimensional variational problem. 

Now that the mathematical problem is formulated, let 
us first have a look at the full solution, which we obtained 
numerically via a conjugate gradient minimization |l2j . 
The results confirm the suspicions made above: For large 
enough aspect ratio £ = ri/ra the equilibrium nematic 
field is n — e„, corresponding to simple circumferential 
winding of the polymer. But as the torus grows suf- 
ficiently fat, a continuous transition occurs to a state 
in which (simultaneously throughout the entire torus) n 
acquires components in e$- and e r -direction, i. e., the 
polymer additionally winds around the tubular circle, see 
Fig. ^ This winding relaxes bending energy, but only at 
the expense of the additional twist, which is zero when 
n = e^. Consequently, this twist instability occurs more 
readily if the ratio r\ = K2/K3 of twist and bend mod- 
ulus is small. All this is confirmed in Fig. |3 where we 
show the maximum twist angle, which is a suitable or- 
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FIG. 2: Order parameter (j> 

'. — &rccos (n^mm) as a func- 
tion of the aspect ratio £ = r\/r2- The curves corre- 
spond to different ratios of elastic moduli, r\ — K2/K3 £ 
{0.01,0.05,0.1,0.2,0.3,0.4}, the gray arrow pointing toward 
increasing values. The inset illustrates the toroidal coordinate 
system. 

der parameter, as a function of the aspect ratio £ of the 
torus. Note that the calculation is reliable down to the 
value £ = 1, where the hole in the torus degenerates to 
a point, the only possible tangent vector is and thus 
4> = it/2. 

The mathematical task of functional minimization can 
often be accurately approximated by devising a varia- 
tional ansatz which is analytically tractable. The follow- 
ing choice turns out to be remarkably good: We will first 
assume that the nematic field does not have a compo- 
nent in e r -direction. The tangential boundary condition 
n r (r — r2) — is then automatically taken care of. The 
remaining two components must satisfy the normaliza- 
tion condition n^ + n 2 ^ = 1, and it thus suffices to specify 
one of them, say n$. It is easy to check that any ansatz 
of the form n$ = /(r)/[l — [r/r\) cost?] with an arbitrary 
function f(r) yields a divergence free nematic field. We 
choose a linear f(r) = Lur/r2, i. e. 

n$(r,d;uj) = u- ^ , (3) 

1 — (r/ri) cosi> 

where ui, which we may call the "twisting strength" , is the 
only free parameter of the ansatz. With this choice we go 
back into the Frank-Oseen free energy calculate the 
derivatives, and integrate over the volume of the torus. 
Since the sign of u> only determines the handedness of 
the twisted structure, it cannot influence the free energy 
E, which thus must be an even function of to |l8l |. An 
expansion in powers of ui 2 then yields 

■J- = go(0+92(^v)^+94(^v)^ + 0(u; 6 ) , (4) 



3 




0.01 - 



FIG. 3: Structural phase diagram of the toroidally wound 
complex on a log- log scale. For small £ = ri/r% and r\ — 
K2/K3 the polymer is wound in a twisted way, for large £ 
and r\ it prefers to wind straight. The dots are results from 
the full numerical minimization, the solid line is Eqn. ©, and 
the dashed line stems from the "improved ansatz" (see text). 



where the expansion coefficients gi are functions of the 
system parameters |14| . In particular, go has the simple 
form 
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This term contributes even if the twist o> vanishes. In 
fact, it coincides with the bending energy of a chain 
of length L which has a curvature energy A/p 2 per 
unit length (p being the local radius of curvature), and 
which is wound without additional looping within a torus 
of volume V — provided the (intuitively clear) relation 
K3V = AL holds. This relates the nematic bending mod- 
ulus K3 to the more usual polymer stiffness A. 

While go helped us to map our parameters, gi will 
localize the phase transition. The reason is that Eqn. Q 
has the form of a Landau free energy as it occurs for 
phase transitions with a scalar order parameter, u>, and 
it predicts a continuous transition at the point where the 
coefficient of the quadratic term vanishes, i. e., 52(^7 v) = 
0. This results in the phase boundary 



V 



e 



1 



4£ 2 



1 + 6 



1 



(6) 



This boundary is shown in Fig. together with exact 
points originating from the full functional minimization. 
It quantifies the conclusion that toroidal condensates will 
spontaneously twist if the aspect ratio £ = r± /r% is small 
(and the torus thus fat), and if the ratio of elastic mod- 
uli T] = K2/K3 is small, i. e. if twisting a bundle is easy 
compared to bending it. For large £ Eqn. © asymp- 
totically behaves like rj ~ Tg£ -2 : thereby explaining the 



exponent —2, which is also indicated in Fig. [3| and which 
actually describes the whole phase boundary quite well. 
The agreement between the simple one-parameter varia- 
tional ansatz and the full calculation is remarkably good. 
It can even be improved by including an additional pref- 
actor 1 — — into Eqn. Particularly for small £ this 
improved ansatz agrees better with the exact answer, ba- 
sically since the new prefactor cancels the unphysical di- 
vergence of the denominator for r — ► n at $ = 0. The 
analytical expression for the phase boundary is quite in- 
volved and will not be shown here, but it has the same 
large £ asymptotics (see Fig. 0). 

Both the ansatz as well as the full numerical solution 
point to an upper critical ratio of elastic moduli, ?/ c , 
beyond which the system will no longer spontaneously 
twist. Even at the lowest possible aspect ratio £ = 1 the 
energetic cost for twisting has become so large that it 
no longer admits a bending relaxation. The ansatz @ 
gives rj c = |, the improved ansatz gives j] c k 0.829, while 
the full numerical solution suggests the deceptively sim- 
ple result r] c = 1. We have no analytical support for the 
latter, but we also want to remind that the limit ^ — ^ 1 is 
somewhat academic, because our tacit assumption that 
the condensate shape is strictly toroidal will most likely 
break down in this case . 

We have thus seen that the ratio r\ and the torus geom- 
etry £ uniquely specify the twist-state of the condensate. 
However, while one can easily measure £ in an experi- 
ment (just by visual inspection), it is hard to specify in 
advance. In contrast to that, the surface tension a can 
be readily changed (for instance via the concentration of 
condensing agents), but its actual value is hard to mea- 
sure. In our simple model it is of course not difficult to 
add a tension term a times the torus surface S to the con- 
densate energy. Using S oc £V3^2/3 ; 

one can re-express 

the twisting transition as being driven by increasing a, 
and it remains continuous. However, practical consider- 
ations would advise to tune a only for the purpose of 
modifying £, but subsequently use £ as the independent 
variable. This way one needs no longer (neither theo- 
retically nor practically) worry about how a particular 
concentration of condensing agents gives rise to a partic- 
ular torus geometry. 

After these local considerations it is time to study 
global aspects of the polymer structure. Let us start 
with the flow itself. It can be shown that incompress- 
ibility, V • n = 0, together with axial symmetry causes 
the flow to be Hamiltonian — hence one more conservation 
law exists . In our case it forces the flow lines to stay 
on two-dimensional slightly deformed toroidal surfaces, 
such that the total flow consists of a nested structure of 
invariant tori. In fact, our ansatz © follows readily from 
the quadratic Hamiltonian H = ^uir 2 , which is constant 
on circular tubular layers. Of course, the actual polymer 
has to switch between these layers, reminding us that ir- 
respective of twist none of the above structures can be 
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realized without localized defects |17| . 

There is one global aspect of the polymer structure 
in which it differs fundamentally from a plainly wound 
torus: As is visible in Fig. 2] the path of the polymer 
threads it repeatedly through the middle hole. Moreover, 
the amount of this looping (as measured e. g. by the av- 
erage change in i? per turn) depends on the layer. This 
effect implies that the entire polymeric strand is heav- 
ily entangled with itself. A rough estimate for £ = 1.5 
and polymer length L = 15/um gives about 30 thread- 
ings through the hole. In fact, were it not for the two 
free ends, these knotted states were topologically inac- 
cessible. In other words, unlike the initial collapse, the 
second stage, the structural ripening, relies on the motion 
of the free chain ends and is thus a much slower process. 

On the other hand this structural ripening meets no ki- 
netic barriers during the relaxation to its twisted ground 
state, downhill the free energy landscape. The motion of 
the two free ends is then energetically directed and their 
local rearrangement does not involve the highly improb- 
able threading through the toroid hole in 3D space. In 
addition, the weak chiral interaction of DNA molecules 
neglected above, gives rise to a (small) chiral term 
in the elastic free energy, which might contribute to the 
symmetry breaking and "guide" the twisting in a pre- 
ferred direction. 

It is also tempting to explain the unusually short DNA 
cholesteric pitch (10 times shorter than in vitro) in sperm 
chromatin assuming that the weak DNA chirality merely 
determines the handedness of twist whereas its pitch is 
given by the twisted state of polymer strands after the 
topological ripening took place. Moreover, the severe 
knotting in highly dense phage heads, where the 2nm 
thick DNA has little space left for usual entropic entan- 
glement effects, suggests the same explanation — in par- 
ticular since it is far more pronounced for genetically 
modified phages in which the DNA is no longer attached 
to the capsid and can thus undergo structural ripening. 

Finally, is the predicted effect strong enough to be 
of some relevance for DNA condensation and gene de- 
livery? For typical experimental parameters of DNA 
length L = 15 — 30/im, £ ~ 1.5 — 2.5, bending stiff- 
ness A ~ 50fceT-nm, and inter-helical distance d ~ 3nm 
we obtain first the elastic constant K3 ~ A/d 2 ~ 20pN, 
which is dominated by the bending stiffness 9]. The 
twist constant Ki can be estimated by the decondensa- 
tion force ~ 2pN obtained in single molecule experiments 
with condenser spermidine [18j. The difference in elas- 
tic energy between the twisted and untwisted states, as 
bounded below by the the variational ansatz, lies in the 
range 15 — SOk^T. This indicates that topological ripen- 
ing stabilizes the condensate. In addition, if the solvent 
quality abruptly improves, the twisted DNA will unfold 
more inertly than its untwisted counterpart, due to heavy 
entanglement with itself. If this stabilization occurs on 



the typical time scales relevant for gene therapeutical ap- 
plications, it might prevent a premature digestion of the 
genetic material by the host organism and influence (pos- 
itively) the efficiency of the gene delivery process. 

We thank H. Schiessel, K. Kremcr, I. Pasichnyk, and 
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